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Abstract 

The Principle of Least Action is used with a simple Lagrangian 
density, involving second-order derivatives of the wave function, to 
obtain the Schrodinger equation. A Hamiltonian density obtained 
from this simple Lagrangian density shows that Hamilton's equations 
also give the Schrodinger equation. This simple Lagrangian density is 
equivalent to a standard Lagrangian density with first-order deriva- 
tives. For a time-independent system the Principle of Least Action re- 
duces to the energy variational principle. For time-dependent systems 
the Principle of Least Action gives time-dependent approximations. 
Using a Hartree product trial wave function for a time-dependent 
many-boson system, we apply the Principle of Least Action to obtain 
the Gross-Pitaevskii equation that describes a Bose-Einstein conden- 
sate. 



1 Introduction 

The Lagrangian approach using Hamilton's Principle of Least (or Stationary) 
Action has been a unifying principle in almost all areas of physics to obtain 
dynamical equations [I]-|3]. The Hamiltonian approach, on the other hand. 
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depends on a Lagrangian to obtain a canonical momentum and Hamilto- 
nian. In nonrelativistic quantum mechanics the Hamiltonian approach has 
dominated, however, because the Schrodinger equation uses the Hamihonian 
operator obtained by quantizing the canonical momentum of the particles. 
The wave function in the Schrodinger equation is a complex function in 
configuration space. From the wave function a Lagrangian density can be 
constructed and used in Hamilton's Principle of Least Action to give the 
Schrodinger equation. This same procedure is used in relativistic quantum 
field theory, which is primarily based on the Lagrangian approach because 
the Lagrangian density is a scalar density [H O El [3] . Our approach is dif- 
ferent from the Feynman path integral formulation of quantum mechanics 
that uses a classical particle Lagrangian in the path integral to obtain the 
propagator [7]. 

In this paper we discuss a simple Lagrangian density for nonrelativistic 
quantum mechanics that gives the Schrodinger equation when used in the 
Principle of Least Action. This Lagrangian density has been called "seem- 
ingly artificial" [S] because it involves second-order spatial derivatives and 
is apparently complex, but it is nevertheless widely used From this 

Lagrangian density we use the canonical procedure for fields to obtain a 
canonical momentum conjugate to the wave function as a generalized co- 
ordinate and then find a Hamiltonian density. Using these, we show that 
Hamilton's equations also give the Schrodinger equation. We then show that 
this simple Lagrangian density is equivalent to a standard Lagrangian den- 
sity that is real and has only first derivatives [3]. In addition, we show 
that for time-independent systems the Principle of Least Action reduces to 
the energy variational principle of nonrelativistic quantum mechanics [T2] . 
For time-dependent systems the Principle of Least Action can be used with 
time-dependent trial wave functions to obtain approximate equations and so- 
lutions. We apply this method to a system of charged bosons with a Hartree 
product time-dependent trial wave function to obtain the Gross-Pitaevskii 
equation [ll]-[l6j, which is a nonlinear Schrodinger equation that describes 
a Bose-Einstein condensate [T7|. 

In Section 2 we review the Hamiltonian for a nonrelativistic quantum 
system of many charged particles in an external electromagnetic field. In 
Section 3 we give a simple Lagrangian density for the system and show that 
Hamilton's Principle of Least Action gives the Schrodinger equation. Using 
this simple Lagrangian density, we apply the canonical formalism in Section 
4 to obtain a Hamiltonian density and show that Hamilton's equations also 
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give the Schrodinger equation. In Section 5 we give an elementary deriva- 
tion of the equation of continuity in configuration space for probabihty based 
on the invariance of the Lagrangian density under infinitesimal global gauge 
transformations. We show in Section 6 that our simple Lagrangian density 
with second-order spatial derivatives is real and equivalent to the standard 
one with first-order derivatives. For a time-independent system we show in 
Section 7 that Hamilton's Principle reduces to the energy variational princi- 
ple. For a time-dependent system we show in Sec. 8 that a time-dependent 
trial wave function may be used in Hamilton's Principle and then apply it 
to obtain a time-dependent approximate equation for a many-boson system. 
The conclusion is given in Section 9. 

2 Hamiltonian for a many-particle system 

For a system of many charged particles in an electromagnetic field the La- 
grangian of the classical system can be used to obtain a particle Hamiltonian 
by using the canonical procedure [IB]. The system is quantized by replac- 
ing canonical momenta conjugate to the coordinates by operators satisfying 
canonical commutation relations [19] and the Hamiltonian then becomes an 
operator. Because of the sum over all particles, the Hamiltonian operator em- 
phasizes the particle aspect of quantum theory. We review the Hamiltonian 
here to establish the background and notation. 

The Hamiltonian operator H for a system of A^-particles with masses 
and charges qi in an external electromagnetic field is 

^ r 1 1 

H = Y.\^. [P-*A(r,, t)f + g.Ao(r„t) \ + V{t^, r^), (1) 

where \^ is a conservative potential energy. The canonical momentum op- 
erator pi= —iKVi = —ihd/dYi is conjugate to the coordinate rj for particle 
i = 1,2, ■ ■ ■, N. Together the coordinates and canonical momenta satisfy 
the canonical commutation relations. The external classical electromagnetic 
field is characterized by a vector potential A{r-,t) with components Aa{r^,t) 
for a = 1,2,3 and a scalar potential Ao{ri,t). From these potentials the 
electric and magnetic fields can be determined in the usual way. 

It is convenient for the rest of the paper to use a more compact notation 
for the Hamiltonian ([T]). The a-component of the displacement of particle 
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i has components Xia for a = 1,2,3. The a-component of the mechanical 
momentum operator Pj for particle i is 

Pia ~ pia Qi^ai (2) 

where the a-component of the canonical momentum operator is pia = —ihd/dxi, 
= —ihdia- Using this notation for mechanical momentum, we can rewrite the 
Hamiltonian operator ([T]) as 

H = ^PL + q.Ao{T,,t) + V, (3) 

where summation over repeated particle numbers i from 1 to iV and repeated 
component indices a from 1 to 3 is understood. 



3 Lagrangian density and Hamilton's Princi- 
ple 

The dynamical equations in almost all areas of physics can be obtained from 
Hamilton's Principle of Least Action [H [2] . For a many-particle system the 
wave function \1/ in nonrelativistic quantum mechanics is a complex func- 
tion in a 3A^-dimensional configuration space. In order to use Hamilton's 
Principle of Least Action to obtain the Schrodinger equation we need a La- 
grangian density. Such a Lagrangian density can be constructed from the 
wave function, its complex conjugate, and their partial derivatives to any 
order [20]. Because it involves the wave function, the Lagrangian density 
emphasizes the wave aspect of quantum theory. The wave function of the 
nonrelativistic quantum system of N particles is a complex, time- dependent 
function \I' = = \l/(r, t) = \l/(ri,r2, ...,rN,t), where the 3A^-dimensional 
vector r = (ri, r2, r^) is a vector in configuration space. 

We postulate a simple Lagrangian density C for this system in terms of 
the nonrelativistic wave function \1/ as 

C = ^*{ihdt- H)-^, (4) 

which depends only on \E'*, \E' and its partial derivatives dt"^ = d'^/dt, dia"^ = 
d^iz/dxia, and df^"^ = d^'^/dx^^. Any Lagrangian density that gives the 
correct equation of motion is a valid one, so a Lagrangian density that is 
complex can still be useful |21j. The Lagrangian density need not have only 
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first-order derivatives as is sometime thought, but can have derivatives of 
any order [20]. Nesbet [S] has called the Lagrangian density (jl]) "seemingly 
artificial," but it is nevertheless widely used [9]-[T3]. The Lagrangian density 
for the Dirac equation is similar in form [1]. 

To obtain the dynamical equations from a Lagrangian density we use 
Hamilton's Principle of Least Action. Hamilton's Principle [3J states that 
the action functional 5'[\E'*, \1'] for any Lagrangian density L is stationary, 

% = j dt j C{^<, ...) = Stationary, (5) 

where the integration is over all time and all configuration space with an 
element of volume dr'^r = d^ri ■ ■ ■ d^r^. The action is stationary when 
its variation with respect to either or \l/ (or both) is zero. For an arbitrary 
Lagrangian density /^(^E^, '^*, ...) variation with respect to gives the Euler- 
Lagrange equation [6]. Boundary conditions at infinity on the variations are 
chosen such that 5^* and 6"^ vanish sufficiently rapid as |t| cxd and rj oo 
for all i. 

When the specific Lagrangian density C in Eq. (jl]) is substituted into Eq. 
and variation is made with respect to \E'*, we obtain 

SS[^*, ^] = j dt j d^^r (5^* (^ihdt -H^^ = 0. (6) 

Since the variation 6"^* is arbitrary except for vanishing at the boundaries, 
Eq. ([6]) immediately gives the Schrodinger equation 

H^> = ^h—, (7) 

where the wave function ^1^(0) must be specified. If the variation of the 
action is made with respect to then integration by parts is necessary and 
the complex conjugate of the Schrodinger equation is obtained. 



4 Hamiltonian density and Hamilton's equa- 
tions 

From the Lagrangian density in Eq. (jl]) we can obtain the Hamiltonian 
density by the canonical procedure. The Hamiltonian density can be used 
to obtain Hamilton's equations. To be consistent with the Principle of Least 
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Action, Hamilton's equations must also give the Schrodinger equation ([7]), 
and indeed they do. This section uses functional derivatives and is more 
detailed for pedagogical purposes. 

The definition of the Hamiltonian from the Lagrangian for classical par- 
ticles may be generalized for fields [18]. The many-particle wave function 
\E'(r, t), where time t is a parameter and r is a vector in the 3A^-dimensional 
configuration space, can be considered as a generalized coordinate with a 
conjugate momentum n(r, t) in the canonical formalism. The Hamiltonian 
density 7i is given by a Legendre transformation [18] to obtain a function of 
and H(t) from the La grangian density C (j4]), 

nm), * W) = n(r, t)^{r, t) - C{^{t), ^{t)), (8) 

where ^l/(r, t) = 9f\l/(r, t). The canonical momentum n(r, t) conjugate to the 
generalized coordinate \l/(r,t) is chosen to eliminate the coefficient of in 
the variation of Eq. ([8]), which gives 

dC 

n(r, t) = = ih^*(r, t). (9) 

Substituting H(r,t) = ih'^*{r,t) and £ in Eq. (g]) into Eq. ([8]), we find 
that the Hamiltonian density is 7i = which is the density of the 

Hamiltonian operator H in Eq. ([T]). The spatial integral of the Hamiltonian 
density is the energy £(t) of the system 

S{t)= j n = j d?^r ^*(r, t)if^(r, t). (10) 

However, when the integral of the Hamiltonian density in Eq. (fTOj) is written 
in terms of the canonical momentum H(r, t) = ih'^*{T, t) it can be considered 
as a Hamiltonian functional EI[\l/(t), H(t)] for the system 

H[^(t), U{t)]= J d^^r n{t) = ^J d^^'r H(r, t)^^(r, t). (11) 

From the Hamiltonian functional and the Euler-Lagrange equation Hamil- 
ton's equations [3l [22] may be written in a familiar form using functional 
derivatives. 

Hamilton's first equation for fields is [5] 

• , , 5e[^, HI , , 
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where the left-hand side is \i/ = dt"^ = d'^/dt. The right-hand side is the 
functional derivative of the Hamiltonian functional EI[\E', 11] with respect 
to the canonical momentum n(x, t),where x = (x-^,X2, ...,XAr) is also a 3A^- 
dimensional vector in configuration space. When the Hamiltonian functional 
in Eq. (|TT|) is substituted into it, Hamilton's first equation gives 

as) 

which is the Schrodinger equation ([7j). The functional derivative of H(r, t) [or 
\E'(r,t)] with respect to H(x, t) [or \E'(x, t)] in Eq. (fT3l) is a 3A^-dimensional 
delta function [6], 

fHI^ -X). (14) 

Integrating over this delta function in Eq. f|T3|) . we obtain the Schrodinger 
equation ([7]). The dimensions of the delta-function on the right-hand side 
of Eq. (JT^ are (length)'"^^ , so the functional derivative also has the 

same dimensions, contrary to appearances. Thus, it is only apparent that 
there is a lack of dimensional consistency in Eqs. fll2l) and (1131) . 
Hamilton's second equation is [5] 

"(^■*)-^' 

where the left-hand side is tl{x, t) = (9jH(x, t) = ihd'^*{-x, t)/dt from Eq. 
The right-hand side is the functional derivative of the Hamiltonian functional 
with respect to the generalized coordinate \E'(x, t). Using the Hermiticity of 
the Hamiltonian, taking the functional derivative inside the integral to obtain 
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a delta function and integrating, we obtain Hamilton's second equation, 



5\E'(x,t) J ih 

(5^(x,t)^ 



^/^^5^^ = -i^^*(x,t), (16) 

which is the complex conjugate of the Schrodinger equation ([7j). There- 
fore, Eqs. (fT3|) and (fT6|) show that Hamilton's equations obtained from the 
canonical formalism using the Hamiltonian density ([8]) obtained from the La- 
grangian density (jlj) is consistent with Hamilton's Principle of Least Action 
in giving the Schrodinger equation. 



5 Conservation of probability 

Noether's theorem states that the invariance of a Lagrangian density under 
a symmetry operation implies a corresponding conservation law [211 [23] • In 
particular the invariance of Lagrangian density (jlj) under infinitesimal global 
gauge transformations gives the equation of continuity for probability [3l [21] . 
Even though this relation is well known, the proof of it at an elementary level 
for nonrelativistic quantum mechanics in configuration space is not readily 
available. We give a short derivation of it here to use in the next section. 
First consider the variation of the Lagrangian density with respect to 
Using "differentiation by parts" [udv = d{uv) —vdu] and the Schrodinger 
equation ([7]), we obtain 

SC = thdt + ^dia UPiaW^^ + ^*PiJ^] ■ (17) 

At this point the variation is arbitrary. When we obtain from an 
infinitesimal global gauge transformation, we find the equation of continuity 
for probability. 

An infinitesimal global gauge transformation on the wave function is 

^' = exp {i6T/h} * = ^ + {i6T/h)^ + ■ ■ •, (18) 

where ST = Sj^j^A is an infinitesimal constant. The vector and scalar po- 
tentials are of course unchanged by a constant gauge function 6A. The new 
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Lagrangian density CJ = exp {—iST / h} C exp {iST / h} = C is obtained using 
the new wave function in place of \l/ in Eq. (jlj), so (5£ = £' — £ = 0. 

The variation in the wave function under the transformation f|T8l) is 
S-qj = -qj' - -qj = [iST/h)-^ [25j. Substituting this 6-^ into Eq. ([T7D and using 
5£ = 0, we obtain a generahzed equation of continuity 

dt [^*^] + ^d,^ [^{P^aW + = 0. (19) 

This equation can be rewritten as the equation of local probability conserva- 
tion in configuration space 

dtp{r,t) + diaJia{r,t) = 0. (20) 

The probability density p(r, t) and probability current density Jia{r, t) in the 
3A^-dimensional configuration space are defined respectively as 

p(r,t) = ^*(r,t)^(r,t), J,,(r, t) = Re ^*(r, t){;,,^(r, t), (21) 

where the velocity operator Vj is the mechanical momentum Pj divided by 
the mass for particle i, i.e., = Pi/rrii. 



6 Standard Lagrangian density 

The "standard" Lagrangian density Ci that gives the Schrodinger equation 
is real and has only first-order derivatives [H |2Tj , 



- -L{p^^q,y{p^^q,) -qf*[V + g.Ao(r,,t)] ^, (22) 

where the subscript on Ci emphasizes that it has only first-order partial 
derivatives. The Lagrangian density Ci in Eq. (122!) can be rewritten using 
differentiation by parts as 

C, = C-z^ U{^*^) + ^d,^ {^*P.a^) I . (23) 



9 



Even though the Lagrangian density £ in Eq. (jl]) is apparently complex, its 
imaginary part is zero, which can be shown by taking the imaginary part of 
Eq. 

Im £ = I !^dt{^*^) + ^dio^Re (^*P.a^)| = 0, (24) 

from the equation of continuity (flQl) . The real part of C in Eq. fl23l) is 

Re £ = £ = - -^9,«Im (^*Aa^), (25) 
znii 

so the Lagrangian densities C and C i are not equal. Nevertheless, they are 
equivalent to each other in the sense that they have the same action [2T] . 
Their difference is a partial derivative with respect to the spatial coordinates 
which vanishes because of boundary conditions when integrated to obtain the 
action The same dynamical equation, viz., the Schrodinger equation, is 
therefore obtained from both Lagrangian densities. 

Sil [26] was the first to use the real part of £ as a Lagrangian density, 

Lsii = Re ^* (ihdt - H^^ = Re C, (26) 

which was later used by Mittleman [27]. It was subsequently used by Mc- 
Carroll, et al. [2S], who pointed out that Re C = C 

7 Energy Variational Principle 

If the quantum mechanical system is time independent then Hamilton's Prin- 
ciple reduces to the energy variational principle of time-independent quan- 
tum mechanics [1^. The time-independent Schrodinger equation can be ob- 
tained by separation of the time variable in the time- dependent Schrodinger 
equation. However, the time-independent Schrodinger equation can also be 
obtained from Hamilton's principle by using it with a trial wave function. 

The Principle of Least Action can be used with a trial wave function of 
the form 

^(ri, rs, ...r;v, t) = ^(ri, rg, ...r^v) exp {-iEt/h - e\t\) , (27) 

where \& on the right-hand side is time independent, E is the energy and 
e > is a small parameter that can be taken to be zero at the end of the 



10 



calculation. Using this wave function in Eq. (jS]) for the action and doing the 
time integration, we obtain 

S[^*, ^] = 1 y ^* (^E - ij) ^ = Stationary. (28) 

If this equation is multiplied by — e < and the constant E is added, we 
obtain the energy variational principle for a stationary state 

j -E{j - 1) = Stationary, (29) 

which is independent of e. The energy ii^ is a Lagrangian multiplier that 
ensures the normalization of the wave function. 

If the variation of Eq. (I29l) is made with respect to the function we 
obtain the time-independent Schrodinger equation 

H"^ = E^. (30) 

The variation with respect to gives the complex conjugate of the time- 
independent Schrodinger equation fl5Ul) after integration by parts. 

Hamilton's Principle of Least action is however more important in obtain- 
ing approximate solutions for the energy and wave function. If the wave func- 
tion \l/ in Eq. (129|) is replaced by a time-independent normalized trial wave 

function $ with parameters or unknown functions and an energy = 

E'. Then Eq. (!29l) can be varied with respect to these parameters or func- 
tions to obtain the lowest energy E' with a wave function of the form 
From the standard Rayleigh-Ritz method it is easy to show that the energy 
^$1^$^ = E' is an upper bound to the true ground state energy Eq, i.e., 
E' > Eq [19]. 

8 Time-dependent approximation for a many- 
particle system 

Approximate solutions or equations for the time-dependent Schrodinger equation 
(I7j) can be obtained from the Principle of Least Action ([5]) by using the La- 
grangian density (jl]) and replacing the exact wave function '^{t) in by a 
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time-dependent trial wave function 

$] = j dt j d^^r(l>*{t) (ihdt - = Stationary, (31) 

and varying with respect to parameters or functions in Equation fl3T!) 
can be used to obtain a time-dependent approximation for any quantum 
system regardless of the number of particles. 

For iV-identical particles with Qi = q and rrii = m, the potential energy V 
in the Hamiltonian ([T]) is taken to be the sum of one- and two-body potentials 

N N 

V{r,,r,,...,r^) = J2v^'\r,) + - V^^'Hr.r,), (32) 

i=l i¥'j=i 

respectively, where the two-body potential is symmetric V^'^^ (rj, r^) = V^'^^ {vj, rj). 
This form of the Principle of Least Action is applicable to any system of many 
identical particles, either bosons or fermions. 

We illustrate Hamilton's Principle with an approximate wave function for 
a system of many bosons to obtain the Gross-Pitaevskii equation [T^-p^. 
which is a nonlinear Schrodinger equation commonly used to treat Bose- 
Einstein condensates [17]. We choose a time-dependent trial wave function 
$(t) for a system of A^-identical bosons as a Hartree product of identical, 
normalized single-particle functions (p{r,t), 

$(ri, r2, rjv, t) = (p{ri, t)ip{r2, t) ■ ■■(/^(rAr, t), (33) 

which satisfies Bose-Einstein statistics because it is symmetric under inter- 
change of particles. Substituting Eq. fl33|) into the variational principle fl3Tl) 
and integrating over the coordinates, we obtain the action functional in terms 
of single boson wave functions 



= N j dt j d^r^*{r,t) !^ihdt^{r,t) - ^ [-ihV-q Af ^{r,t) ~ q Ao^{r,t) 

-V«(r)¥,(r,t)-^(iV-l) I dV02)(r,r')¥.*(r',t)^(r',t)y,(r,t)| 

= Stationary, (34) 
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where the integration is over all time and space. 

Variation of this expression with respect to (p* gives a time- dependent 
Hartree approximation known as the Gross-Pitaevskii equation [H]-|16j. 

z;i^(/.(r, t) = ^ l-thV-q A(r, t)f ^{r, t) + [q Ao{r,t) + V^^\r)]ip{r, t) 
+ {N-l) j rfVr(2)(r, v')ip*{v\ t)ip{v', tMr, t), (35) 

for a charged many-boson system in an electromagnetic field [29]. This equa- 
tion is a nonlinear Schrodinger equation for the single boson wave function 
ip, where the nonlinear term is the average potential on one particle due to 
all the other particles in the system. The variation may also be made with 
respect to ip{r,t) to obtain the complex conjugate of Eq. fl35|) . 



9 Conclusion 

Hamilton's Principle of Least Action is a fundamental principle of physics 
that is used to obtain dynamical equations for both nonrelativistic and rela- 
tivistic particles and fields [U E] . This principle is applied to nonrelativistic 
quantum mechanics by considering the wave function in configuration space 
as a generalized coordinate and constructing a Lagrangian density such that 
Hamilton's Principle gives the Schrodinger equation. This approach gives 
a unified treatment for both nonrelativistic quantum theory and relativistic 
quantum field theory [1], as well as showing their unity with other branches 
of physics. The Lagrangian approach is in a sense more fundamental than the 
Hamiltonian approach that has dominated the treatment of nonrelativistic 
quantum mechanics because the Hamiltonian must be derived from a La- 
grangian. In relativistic theory the Lagrangian density is also a scalar density, 
whereas the Hamiltonian density is not. The two approaches are, however, 
complementary in the sense that the Lagrangian density using the wave func- 
tion emphasizes the wave aspect of nature, whereas the Hamiltonian operator 
with a sum over particles emphasizes its particle aspect. 

We use a simple Lagrangian density for the Schrodinger equation that 
is apparently complex and has second-order spatial derivatives. This La- 
grangian density is actually real because of the equation of continuity for 
probability. It is equivalent to the standard Lagrangian density that is real 



13 



and has only first-order derivatives [21] . Two Lagrangian densities are equiv- 
alent if the difference between the two is a sum of spatial or temporal partial 
derivatives so they give the same action. 

For time-independent quantum systems, the Principle of Least Action 
reduces to the energy variational principle of nonrelativistic quantum me- 
chanics. By minimizing the energy calculated from a trial wave function, it 
gives an approximate wave function and an energy that is an upper bound 
to the exact ground state energy. 

For time- dependent quantum systems, the Principle of Least Action with 
a time-dependent trial wave function gives an approximate equation for the 
trial wave function. We apply this approach to a many boson system by 
using a Hartree product wave function to obtain the time-dependent Gross- 
Pitaevskii equation [T1]-[T6], which is a nonlinear Schrodinger equation that 
describes a Bose-Einstein condensate |T7]. For fermion systems a Slater de- 
terminant with time- dependent single-particle trial wave functions may be 
used to obtain the time-dependent Hartree- Fock equations [8]. 



The Lagrangian approach to nonrelativistic quantum mechanics described 
here would be suitable to supplement a course in quantum mechanics or field 
theory. For quantum mechanics it would show how to obtain time-dependent 
approximations. For field theory it would be a simple familiar example. For 
both its use would show the unity that Hamilton's Principle of Least Action 
provides for physics PP-[3]- 
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